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Abstract. The incipient quantum phase transitions of relevance to non-zero 
fluctuations and entanglement are studied in Heisenberg clusters by exploiting the 
negativity as a measure in bipartite and frustrated spin-1 anisotropic Heisenberg 
clusters with bilinear-biquadratic exchange, single-ion anisotropy and magnetic field. 
Using exact diagonalization technique it is shown that quantum critical points signaled 
by qualitative changes in behavior of magnetization and particle number, is ultimately 
related to microscopic entanglement and collective excitations. The plateaus and 
peaks in spin and particle susceptibilities, define the conditions for a high/low density 
quantum entanglement and various ordered phases with different spin (particle) 
concentrations. 
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1. Introduction 

Entanglement properties for few spins or electrons can display the general features of 
large thermodynamic systems and different measures of entanglement have been defined 
to understand QPTs [U E] . A finite spin system is important in the context of molecular 
magnetism and spin pairing. A new line of research points to a connection between the 
local entanglement in one-dimensional correlated systems and the existence of a QPTs 
and scaling [U |5] relevant to the quantum critical points (QCPs). Furthermore, such 
a connection can be exploited to unveil a fundamental connection between the QCPs 
in finite-size small, large clusters (6j [7] and macroscopic systems A particle and 
spin density fluctuations, extending the essential properties of entanglement beyond the 
conventional framework, have been introduced with the explicit reference to the phase 
transitions in canonical and grand canonical ensembles signaled by a critical behavior 
in terms of the energy gaps and susceptibilities [H |9]. The quantum gas of clusters 
at the equilibrium gives unprecedented opportunity to explore exactly these ideas for a 
quantum dynamics of spin fluctuations [10] • While the basic features of entanglement in 
spin-i systems are by now fairly well understood [H], entanglement properties of larger 
spin fermions (or bosons) are less known due to the lack of good operational measures for 
high spin entanglement [12]. A general classical spin-1 Blume- Emery- Griffiths (BEG) 
model [33] has proven to be a useful for description of liquid-gas, liquid-crystal phase 
transitions, tricritical and A points, spontaneous phase separation [HJ [151 HSJ E] • The 
integer spin Heisenberg model exhibits a characteristic spin gap and very rich phase 
diagrams [T8| [T9j . Exact calculations of thermodynamic and entanglement properties 
in finite-size clusters can give an appealing alternative to get insight into the general 
features of bipartite and frustrated systems [20]. Some analytical and numerical studies 
of entanglement and negativity in bilinear-biquadratic spin-1 Heisenberg model on 
dimerized bipartite and frustrated systems have been performed in j22j |23l GU 125] . 
The entanglement with bilinear-biquadratic Hamiltonian has been considered for the 
case of two spins (qubit) in the absence of crystal field [26]. One of interesting problems 
concerning entanglement is to study the effect of uniaxial single-ion anisotropy and 
magnetic fields on negativity. An exact calculation of entanglement versus longitudinal 
crystal field and biquadratic coupling for analyzing the variation of negativity versus 
parameters of the spin-1 system have not been attempted either for ferromagnetic or 
antiferromagnetic Heisenberg model even for small bipartite and frustrated clusters. The 
aim of this work is to discuss, in a general framework, how microscopic entanglement in 
the two- and three-qubit context can be related to a QCPs characterized by plateaus in 
peak behavior of the spin revealed in saddle point singularities on model parameters. We 
provide a reinterpretation of the spin and particle susceptibilities near quantum critical 
points in terms of the quantum entanglement in a physically transparent way. Here we 
adopt negativity to measure the ground state entanglement for spin-1 systems, to reveal 
QPTs in terms of negativity. We have two main goals in this paper: The first is to 
provide a global view of the most general spin-1 Heisenberg model which have not been 
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highlighted so far in minimal clusters. The second is to show that quantum entanglement 
exhibits the existence of characteristic plateaus in negativity related to QPTs. In this 
paper we perform exact calculations of entanglement and response functions in spin-1 
Heisenberg model with bilinear-biquadratic exchange interactions in longitudinal crystal 
and magnetic fields. The basic principles for calculation of negativity are introduced in 
Sect. |3l The ground state magnetic and entanglement properties in spin-1 Heisenberg 
model for ferromagnetic and antiferromagnetic couplings are given in Sect. 14.11 The 
negativity analyzes in absence and presence of magnetic field are given in Sect. 14.21 and 
Sect. |4]3] respectively. The effect of nonlinear interaction is studied in Sect. 14.4]) . Results 
for frustrated trimer are presented in Sect. 14.51 The conclusion is given in Sect. [5] 



2. Spin-1 Heisenberg model 



We consider the spin-1 isotropic Heisenberg model in the presence of magnetic field 
B < 



N 



H = J2[J(SiS i+ i) + K(SiS i+1 ) 2 } + 



i=i 



JV 



j=l i=l 

The linear J and nonlinear K terms are the exchange and quadrupolar interactions. 
Here we implemented the longitudinal crystal field D, which describes an uniaxial single- 
ion anisotropy. In what follows, the crystal field significantly changes the results on the 
entanglement. Notice, effective spin Hamiltonian fl2]) can be derived from Bose-Hubbard 
model in the strong coupling limit. The local spin vector Si for each site has components 
of the Spin-1 operators 
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Unless otherwise specified, we will consider periodic boundary conditions, such that 
Stv+i = Si, where N is the total number of lattice sites. The sum over lattice sites for 
crystal field term with (S 4 2 ) 2 in (j2j) can be reduced to the spin concentration (particle 
number) , 



JV 
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where Pq the number of lattice sites with S- = 0. Notice, the axial anisotropy in many 
respects is analogous to the chemical potential D = —ji. 
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3. Definitions and Basic 

At thermal equilibrium, the state of the system is determined by the density matrix 

H_ gj 

P k B T » — P k B T 

pCO = ^ = ^^H^I, (3) 



where E\ are the eigenvalue of the i-th quantum many body eigenstate and the partition 
function is Z = Yli e ~ l3Ei with = l/fc^T = 1). The many-body entanglement is 
described by the density operator in [201 E3 EHl EH]. For spin-1 system the degree of 
pairwise entanglement, measured in terms of the negativity iVe, can be employed to 
evaluate the thermal state of concern [30J. The negativity of a state p is defined as 

Nc = Y^\i*\, (4) 

i 

where /ij s are negative eigenvalues of p Tl and T\ denotes the the partial pairwise 
transpose with respect to the first system, i.e., for bipartite system in state p it is 
defined as 

(h,j2\p Tl \hJ2) = {kij2\p\ii,k), (5) 
for any orthonormal but fixed basis. Definition (Tj0) is equivalent to 

Ne=M^l, (6) 

where ||p Tl ||i is trace norm of p (p = Tr^/p^p). For unentangled states negativity van- 
ishes, while iVe>0 gives a computable measure of thermal entanglement. 
As thermodynamical characterisation we have used the responses of the thermodynam- 
ical potential with respect to D and B which are follows 

r)F r)F 

P = ^ sz ) 2 ) = w {sz) = dB- (7) 

Here F is the free energy F = —TlnZ and (....) indicates averaging performed within 
a canonical ensemble. The responses for the first derivatives of the thermodynamic 
potential with respect to D and B provide exact expressions for particle xd and spin 
Xb susceptibilities: 

dP d (S z ) 



4. Results 



Phase transitions and entanglement properties in spin-1 Heisenberg clusters with single ion anisotropyh 



4-1. Entanglement and magnetic properties of spin-1 isotropic Heisenberg dimer 

In this section, we consider Hamiltonian in case of N = 2, namely Heisenberg model. 
In the two-qubit case, we diagonalize the Hamiltonian, and obtain the eigenvalues 

E 1 = - 2(B - J - K - D), E 2 = -2(J - K - D), 

E 3 = 2(B + J + K + D), E A = -B - 2 J + 2K + D, 

E 5 = B - 2 J + 2K + D, E 6 = -B + 2J + 2K + D, 

E 7 = B + 2 J + 2K + D, E 8 = - J + 5K + D - A 

E 9 = -J + 5K + D + A , (9) 



and corresponding eigenvectors 
l^i) 
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(10) 
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According to 



2{J-K) 

-1,0,1 and j = —1,0,1) are the eigenvectors of SfSf +1 
Schmidt theorem \ip 5 ) and are not entangled and the maximum entangled states 
can be only \ip 8 ) or \ip$). The partial transpose density matrix of the thermal state p(T) 
at equilibrium is 
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where 
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, 1 B-2{J+K)-D / 4J- 

f — -e t 1 ± e t 
s 2 V 

, 1 2{J-D-K) 

= ± 2 e ^~ + 

e I^R ^ cogh *l + (J_ K _ D ^ sinh Ao 



SI 



A 



2A 

>v, _ 

T 



T_ 



2e d ^r^{K-J) sinh^ 
Ao ' 
e J ~ 5 r~ D (A cosh ^ - ( J - K - D) sinh ^ 



Ao 

here the partition function is 

2{D+J) + 5K D+3K , 4J R 

Z = e T (2e t (l + e T )cosh(-) + 

D+3J /Aoxx 

e T +2eT cosh(— ) + 2e T cosh(— )). 
4-2. Spin-1 in zero magnetic field 

Here we analyze the effects of crystal field on the ground state entanglement in the 
spin-1 Heisenberg model (j2J) at zero field (B = 0). The (quadrupole) particle number 
and negativity plots in figures Q] a and b are both asymmetric as function of D for 
ferro J > and antiferromagnetic J < couplings. The monotonic behavior of P 
versus D in figure Q] a signals a smooth character of the phase transition. Note that 
this smooth bosonic behavior of spin concentration P versus D is contrasted from the 
(step-like) abrupt fermionic change in the electron number as a function of the chemical 
potential in [7j. At infinitesimal T — > the variation of negativity figure db versus D 
for antiferromagnetic case (J > 0) is non monotonic. So, for D = 0, we have a highest 
possible entanglement and ^8 is a ground state of the system. The system for J < 
displays two distinct phases: one separable and other entangled. For positive D region 
the negativity for J > is more than for J < 0. For D — the ground state is five 
fold degenerate (i.e., it is a mixture of ipi,^,^,^,^ states) with zero negativity. At 
infinitesimal D — > +0 the system is entangled in the pure state, For D < the 
ground state at J < is double degenerate with a mixture of V>i and ip3 states. 

(a) (b) 



Figure 1. The density variation of a) the particle number P and b) 
the negativity Ne versus D for antiferromagnetic, J = 1 (dashed) and 
ferromagnetic, J = — 1 (solid) cases 
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Notice, these states are separable (can be factorized) , and therefore, according to 
definition, these quantum states are without entanglement. Thus, the entanglement 
in D < region can be used to detect quantum correlations in antiferromagnetic case, 
which are absent for "classical" ferromagnetic. 




Figure 2. The density of negativity Ne versus of J and D. The crystal field D 
enhances the entanglement at J < 0. 

The negativity is non monotonous function with one maxima at D = for J > 
and in a close vicinity to origin at J < 0. The magnetic and quadrupole susceptibilities, 
i- e Xb Xd, allows to distinguish the ordered and disordered phases in the case of broken- 
symmetry at QPTs. Figure [2] shows the pure (extremal) and mixed (non-extremal) 
quantum states. Disentangled dark region in ferromagnetic case corresponds to the 
plateau-like behavior in zero (spin) magnetic susceptibility xo — ^~|b^o versus J and 
D plane in figure [3] a at J < 0. 

(a) (b) 



j j 

Figure 3. The densities for a) zero field (magnetic) spin susceptibility and b) particle 
susceptibility versus J and D. 

The high density magnetic (spin) susceptibility in white sector corresponds to the low 
density of negativity in figure [2j The strong enhancement of negativity along the line 

D = is relevant to the observed peaks in the particle susceptibility, xd = 9D in 
figure [3] a. The various regions seen for (density) negativity in figure [2] are reproduced in 
density of quadrupole susceptibility in figure [3] b versus D and J. Similarly, the phase 
diagram in K—J space in the absence of D and B fields shows the degree of entanglement 
and phases due to effect of nonlinearity on the eigenvalues and eigenvectors in (TlOl) . For 
example, the J = K line separates the maximum entangled and non-entangled phases 
for ferromagnetic coupling, while the J = 3K line is a boarder between entangled and 
new less-entangled phases for antiferromagnetic coupling. 
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(a) (b) 




B 



Figure 4. The density plot for negativity dependencies on crystal D and 
magnetic B fields at T = for (a) J = —1 and (b) J = 1. For both 
(antiferromagnetic and ferromagnetic) cases there are more than three phases, 
which indicates the possible existence of triple or tricritical points. 

We find that for K > 0, the line J = as before separates non entangled and maximum 
entangled phases. The maximum entanglement, which exists for J < at K < J 
and for J > at J > 3K, corresponds to observed condition for bose condensation of 
unpolarized Na atoms on optical lattice [2T] . 

4-3. Effects of magnetic field 

Magnetic field B partially removes the ground state degeneracy and in figure H] one can 
see the presence of new phase boundaries. The entanglement properties of the excited 
states are independent from those in the ground states. Also we found that the pairwise 
entanglement decreases from ground state to excited states, i.e., the more excited the 
system, the less the entanglement. In ferromagnetic case for D = and B = point 
there is a maximum entangled state. In figures H] a for J = — 1 and H]b for J = 1, the 
energies are measured with respect to \J\, which is set to 1. When D <\B\ the system 
is in ipi or ipz state, which is non-entangled. For fixed B the two consecutive phase 
transitions take place by increasing D at D — \B\ and D = \J\ + 6|-B| + B 2 — 1, into tpQ 
and ground states correspondingly. For antiferromagnetic case, the phase diagram is 
more complex. The negativity contains the triple point at \B\ = | and D = — |, which 
implies the presence of various phases, possible coexistence or phase separation in spin-1 
system. When D < — |, the line \B\ = — ^ 2 +\b\ — ^ se P ar& tes and i-e. maximum 
entangled phase from non-entangled one. For D > — | there are three phases in the 
ground state: non-entangled state at D < \B\ — 4; the maximum entangled between 
l+y/B 2 + 2\B\ -7 and 1- ^B 2 + 2\B\ - 7; non saturated entanglement for ip^ states. 
In figure H] there is no any entanglement beyond some critical field B c restricting the 
black region. Also note that entanglement increases with D. Positive D values favor 
to the larger entanglement, while D < shows the tendency toward non-entangled 
states with larger total spin. The ground state diagrams figures [2, and H] exhibit 
quantum critical behavior on the borderlines between various states with continuously- 
varying quantum critical points separating antiferromagnetically ordered distinct phases 
from from the non-entangled state (spin liquid phase). These critical lines, similar to 
continues QCPs, can be used for the classification of the many-body ground states 
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of interacting spins and quadrupole moment in multidimensional parameter space. 
Dynamic interactions between the spins strongly renormalize various parameters in the 
effective Hamiltonian and, therefore, spin and quarupole momentum have properties 
different from a quasiparticle description. As in QCP [2], various states along quantum 
critical boundaries here are necessarily separated by second order phase transitions 
with various entanglement and susceptibility. The quantum critical (lines) boundaries 
appear to be a useful for understanding the formation of various thermodynamic phases 
in the ground state. These continues boundaries in thermodynamic phase diagram 
at infinitesimal temperature coincide with the corresponding QCPs, derived from the 
peaks of magnetic susceptibilities in agreement with our preliminary analysis (see also 
[HI [8]). The boundaries between the various phases are useful for understanding also the 
behavior of thermal negativity for departure to non-zero temperatures. 

(a) (b) 



(c) 



(d) 




Figure 5. The density plots for (a) particle number P (b) (quadrupole) particle 
susceptibility (c) magnetic susceptibility and (d) negativity versus B and D 
when K = 2 in antiferromagnetic case J = 1. 

The distances along the magnetic field in figures H] between various phases define the 
stable magnetic phases with distinct spin gaps configurations, characterized by different 
spin concentration and diverging susceptibilities along the boundaries. For example, the 
negativity for white areas reaches the maximum (saturated) value, while there are also 
different distinct areas with partial (unsaturated) entanglement. 

(a) (b) 




Figure 6. The density variation of a) particle number P and b) negativity Ne 
versus D for K = (solid) and K = 2 (dashed) for antiferromagnetic (J = 1) 
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in zero field B = 0. 

These density plots can be used to determine of QCPs and the boundaries for various 
QPTs. This result for finite size clusters can have important consequences in the physics 
of quantum phase transitions [2], where so far the usual method to detect a phase 
transition is to look at the scaling in the thermodynamic systems. The competition 
among the different phases can lead to complex behavior with the two triple points. 
The density of negativity is an efficient indicator of QPTs. In figure IHa) we find the 
new spin phase boundaries from black to grey region with jump ^ and from grey region 
into the two white ones with the same jump. The white middle line B = in figure [5](c) 
corresponds to classical effect at J < case (without change in entanglement). On the 
other hand, the continuous lines seen in the same plot at J > correspond to quantum 
phase transitions (observable also in negativity). 

4-4- Effect of quadruple term 

Here we display effect of nonlinear interactions between the spins. The variation P 
versus D is shown in figure |6K for two quadruple interactions K > for antiferromagnetic 
case with J = 1. At K = 2, an opposite spin pairing gap is opened at P = 1/2. 
Such a density profile, showing finite leap near P = 1/2, resembles the MH plateau 
behavior for the number of particles versus chemical potential in the Hubbard clusters. 
This is indicative of a possible opposite spin pairing instability [7J. Therefore, the 
cluster behaves at large K as a MH like insulator in contrast to the spin liquid like 
behavior with the zero gap, shown at K = in figure [6b- As it is seen from the density 
plot the magnetic field and quadrupole interaction (K) makes the phase structure in 
antiferromagnetic case more richer. Our analysis shows that the negativity in D — K 
space for ferromagnetic coupling is always less than for antiferromagnetic case. 




Figure 7. Density plot of negativity via J and D at K = B = in three-sites 
cluster. 

4-5. Non-bipartite clusters 

Finally, we display in this session the results of negativity versus D and J for frustrated 
three-site clusters in figure [7] at rather low T = 0.01. This picture for ferromagnetic case 
resembles corresponding figure [2j However, there is an apparent difference in behavior 
for the region J > 0, where there are two extra continuous boarder lines. 
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5. Conclusion 

In this paper we adopted the concept of entanglement to analyze behavior of the 
small-size spin-1 Heisenberg clusters. We used a negativity as a computable measure of 
entanglement to perform extensive calculations of the negativity and response functions. 
The critical fields and intrinsic parameters beyond which entanglement disappears are 
calculated. We found regions where the quantum entanglement can be increased more 
rapidly by increasing both, D and K. The negativity can determine the boarders 
between ordered phases with excess correlations, above the classical ones. The observed 
plateaus and peaks in spin and particle behavior and susceptibilities can be considered 
as a possible universal method for the simultaneous detection of quantum and classical 
phase transitions. The (density) plots are convenient (topographic map) tool for 
observation of the quantum phases and quantum transitions. The states with vanishing 
classical correlations but existing quantum correlations in entanglement open up the 
new opportunities of phase transitions that are detectable only through correlation in 
behavior of entanglement and thermodynamic properties. Our studies of QCPs in small 
size spin clusters appear to be generic to large thermodynamic systems. The exact 
diagonalization is completely unbiased for the study of QPTs and QCPs in strongly 
correlated spin and electron systems [8]. Although the exact studies have limitations 
(since the computations grow exponentially with cluster size), we do not find a minimal 
critical length in clusters below which a quantum critical behavior disappears. The 
spin-1 boson Hubbard like model at certain conditions can be mapped onto the spin-1 
Heisenberg model. Then these studies can also be useful for the analysis of spontaneous 
phase separation and the transition from Mott insulator to quantum superfluid in spin-1 
BoseHubbard models on optical lattices. 
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